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ABSTRACT 



Using Chandra imaging spectroscopy and VLA L-band maps, we have identified radio galaxies at 
P(iaghz) > 3 x 10 23 W Hz -1 and X-ray point sources (XPSs) at h(o.3-SkeV) > 10 42 er g s s ± m 
11 moderate redshift (0.2 < z < 0.4) clusters of galaxies. Each cluster is uniquely chosen to have 
a total mass similar to predicted progenitors of the present-day Coma Cluster. Within a projected 
radius of 1 Mpc we detect 20 radio galaxies and 8 XPSs (3 sources are detected in both X-ray and 
radio) confirmed to be cluster members above these limits. 75% of these are detected within 500 kpc 
(projected) of the cluster center. This result is inconsistent with a random selection from bright, red 
sequence ellipticals at the > 99.999% level. We suggest that these AGN are tr iggered somehow b y the 
intra-cluster medium (ICM), perhaps similar to the Bondi accretion model of lAllen et al.l (|2006j ). All 
but one of the XPSs are hosted by luminous ellipticals which otherwise show no other evidence for 
AGN activity. These objects are unlikely to be highly obscured AGN since there is no evidence for 
large amounts of X-ray or optical absorption. One XPS, in addition to possessing a pure absorption- 
line optical spectrum, has a large excess of light blueward of the Ca II H&K break that could be 
non-thermal emission; a second XPS host galaxy probably has excess blue light. The most viable 
model for these sources are low luminosity BL Lac Objects, similar to the high-energy-peaked BL 
Lacs (HBLs) discovered in abundance in serendipitous X-ray surveys. The expected numbers of lower 
luminosity FR 1 radio galaxies and HBLs in our sample converge to suggest that very deep radio 
and X-ray images of rich clusters will detect AGN (either X-ray or radio emitting or both) in a large 
fraction of bright elliptical galaxies in the inner 500 kpc. Because both the radio galaxies and the 
XPSs possess relativistic jets, they (and, by extensio n, the entire RLF) c an inject heat into the ICM. 
Using the most recent scalings of Pj e t oc L° 5 from iBirzan et all (|200S), radio sources weaker than 
our luminosity limit probably contribute the majority of the heat to the ICM. Also, because these 
heat sources move around the cluster, AGN heating is distributed rather evenly. If a majority of ICM 
heating is due to large numbers of low power radio sources, triggered into activity by the increasing 
ICM density as they move inward, this may be the feedback mechanism necessary to stabilize cooling 
in cluster cores. 

Subject headings: galaxies: active - galaxies: clusters: general - radio continuum: galaxies - X-rays: 
galaxies: clusters - X-rays: galaxies 



1. INTRODUCTION 

Observations of galaxy clusters have provided impor- 
tant clues about cosmology, structure formation and 
the evolution of galaxies. As the largest gravitationally 
bound objects in the Universe, clusters are unique loca- 
tions to study the creation and evolution of AGN. Based 
upon the observed correlation betwe en supermassive 
Black Hole mass and galaxy bulge mass (jMagorria n et al.l 
1998), we expect that there should be numerous, lumi- 
nous AGN in the bright ellipticals in rich clusters which 
can affect their surroundings. Indeed, observational 
evidence supports the notion that non-gravitational 
processes have affected the entropy of the intraclus- 
tcr medium (ICM), as suggested by the steepness of 
the L x -T x rel ationship compared to self-similar scaling 
expectations (lEdge fc StewartJ 119911 : iMarkevitchl 119981 : 
Ruszkows ki et al.ll2004|) . Radio- loud AGN have been di- 
rectly implicated in cluster ICM heating because evacu- 
ated "bubbles" in the diffuse X-ray emitting gas have 
been observed that are spatially coincident with non- 
thermal radio emission (jFabian et al.ll2002t IB irzan et alJ 
2004). However, up until now only AGN in the bright- 



est cluster galaxies (BCGs) have been observed to spawn 
these bubbles and so only the central AG N has been in- 
corpo rated into ICM heating models (e.g.. lBruggen et"aTI 
120051 ). In this case, theoretical difficulties arise when at- 
tempting to distribute t his heating throughout the cen- 
tral cluster regions (e.g.. IVernaleo fc Reyn olds 2003). A 
natural solution would be to have other cluster AGN in- 
jecting energy thro ughout the inner cluster region 
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Most cluster AGN studies have ei ther employed clus- 
ters in a flux-limited survey (e.g., IStocke et alJ 119991 : 
iRuderman fc Ebelingi 120051 : iBranchesi et all 120061 ) or 
simply selected clu sters based upon availability (e.g., 
iMartini et al.ll2007Q . But since flux- limited surveys pro- 
duce cluster samples which contain more luminous and 
thus more massive objects at higher redshifts, this selec- 
tion naturally identifies more evolved structures at ear- 
lier times. This is clearly opposite to the type of selection 
one would prefer to use to study cluster, galaxy or AGN 
evolution. Moreover, this mismatch can result in com- 
parisons that can obscure any true evolution, since more 
massive clusters are placed at the beginning of the evo- 
lutionary sequence and less massive objects at the end. 
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In this paper, we choose galaxy clusters by using a 
new selection method that avoids the difficulty just de- 
scribed. In this summary of first results, we use this sam- 
ple to investigate the nature of radio galaxies and X-ray 
point sources (XPSs) in clusters at moderate- z (0.2-0.4). 
For now, we avoid using the term X-ray AGN until the 
AGN nature of the XPSs is proven (or disproven) by our 
observations presented below. Previous studies simply 
assumed that these X-ray e mitters are AGN (e. g ., the 
"optically-dull X-ray AGN"; IMartini et aD [200l . 120061 : 
lEastman et~aT1 120071 ) . The AGN nature of these XPSs 
has not been carefully scrutinized before and their X-ray 
and radio properties could have a significant impact on 
models for AGN heating in clusters. We address each 
of these points in detail in this paper. We describe our 
sample selection in § [2] and data analysis of X-ray and 
radio selected AGN in § [31 In § 2] we discuss the obser- 
vational properties of the radio and X-ray cluster source 
populations, including location within the cluster and X- 
ray-to-radio flux ratios. In §[5] we discuss the nature of 
the cluster XPSs. In § [5] we discuss the implications of 
this study, particularly for the cluster heating problem. 
In § we summarize our results and suggest testable 
predictions. We use Ho = 70 km/sec, = 0.70, and 
= 0.3 throughout. While we present a summary of 
our general methodologies and detected sources here, a 
m ore detailed discus sion will be presented at a later time 
in lHart et all (|2009t) . 

2. THE "ROAD TO COMA" SAMPLE 

Guided by hydrodynamical simulations to track the 
growth of massive cluster halos (M > 1O 15 M ) from 
z ~ 1 to the present epoch, we have used the tempera- 
ture of the ICM as a proxy for cluster mass. We selected 
0.2 < z < 0.4 clusters that are predicted to evolve into 
objects like the presen t-day Coma Cluster (kT=8.2±0.2; 
iLutovinov et al.l 120081 ). The details of the cosmological 
simulations (similar to lEvrardlll988D w hich include pre- 
heating of the ICM (iBialek et al.ll200l[) will be fully de- 
scribed in lHart et al.l ( 2009T T In short, measurements 
of Tx(z)/Tx(z=0) for individual massive cluster halos 
in the simulation are normalized by the Coma Cluster's 
present temperature to illustrate the variance of cluster 
properties. Figure [I] displays the median temperature 
of these simulated Coma Cluster progenitors as a func- 
tion of redshift (solid line) and the 25th and 75th per- 
centiles about this distribution (dashed lines). At z=0.3, 
the ICM temperature of a Coma Cluster progenitor is 
predicted to be 7.0±2.6 keV. The temperature spread 
is due partly to the modest number of realizations of 
Coma Cluster mass-scale objects in these simulations. 
At z=1.0, the permitted temperature range for our sam- 
ple slowly decreases to 4.0±2.0 keV. Our unique method 
attempts to select clusters of similar mass at a given red- 
shift and thus these clusters on the "Road to Coma" will 
be used for consistent evolutionary comparisons later. 

We selected clusters readily available in the Chandra 
archives with adequate exposure time to detect X-ray 
sources with L .3-s.okeV ~ 1 x 10 42 ergs s -1 . This lu- 
minosity is well above the expected emission from low- 
mass X-ray binaries in ellipticals (Lx < 5 x 10 38 ergs s _1 ; 
iKim fc Fabbianoll2004D and high -mass X-r ay binaries in 
starburst galaxies (Lx < 5 x 10 40 ergs s~ 1 ; lGrimm et al.l 
120031) . Also, for this redshift range, a 1 Mpc radius region 



falls on one ACIS chip and so has been selected consis- 
tently as our survey radius centered on the peak of the 
diffuse X-ray emission. We also require radio observa- 
tions to detect radio galaxies with Piaghz > 3 x 10 23 
W Hz -1 across the entire survey. Along with our well- 
defined cluster sample, our multi-wavelength approach to 
identi fy cluster AGN differs from many previous studies 
fe.g., iLedlow fc Owen|[l995l IMartini et all 120071) which 
use X-ray or radio observations, but not both, to de- 
tect these objects in only one waveband. Table [T] lists 
our eleven z=0.2-0.4 "Road to Coma" clusters, their red- 
shifts, Chandra observational details (ObsID, ACIS aim- 
point, exposure time, ICM temperature, flux and lumi- 
nosity limits), as well as VLA 1.4 GHz radio details. The 
following section details the X-ray, radio and optical ob- 
servations and data analysis. 

3. MULTI- WAVELENGTH DATA ANALYSIS 

3.1. X-ray Imaging Spectroscopy 

We re-processed the Chandra archival observations us- 
ing CIAO v3.3 and CALDB v2.2. Following the typical 
ACIS data preparation pipeline, event-1 files were repro- 
cessed to remove bad pixels, afterglow pixels, and streak- 
ing patterns. The newest calibration files for charge 
transfer inefficiency and time-dependent gain were ap- 
plied, and then finally the event files were filtered on 
status and grade. Flaring events were identified by ex- 
tracting the 0.3-12 keV count rate on either the ACIS-S1 
or ACIS-I1 chips and we eliminated time periods with 
rates >3cr above the observed mean. 

3.1.1. Cluster ICM Temperature 

Cluster ICM spectra were extracted from within a 1 
Mpc projected radius after obvious point sources were 
removed (point source identification is described below) . 
If a background area could not be extracted on the same 
chip as the cluster emission, we extracted background 
spectra from re-projected blank sky observations. Obvi- 
ous cool-core regions were excised prior to spectral mod- 
eling. Source spectra were binned to a minimum of 20 
cts/bin, then modeled using XSPEC (fA rnau dl fl99fih and 
MEKAL models with foreground extinction held fixed at 
the Galactic hydrogen col umn at the cluster coordinates 
(jDickev fc Lockmanlll990f ) and a constant 0.3 Solar metal 
abundance. 

The lower redshift limit of our cluster sample is lim- 
ited by the size of one ACIS chip, allowing a ~ 1 Mpc 
radius to fall on one chip. A 1 Mpc radius corre- 
sponds to ~ 0.5 R2oo- However, our simulated cluster 
temperatures were extracted from within R2oo- Based 
on the similarity in the temperature profiles of many 
clusters for r>0.1 R2nn, in both cool-core and non-cool 
cores (e.g.. iDe Grandi fc Molendl l2002t IVikhlinin et all 
120051: iBaldi et al.ll2007t iPratt et al.ll2007l ). we expect the 
ICM temperature profile to be decreasing beyond 0.5 
R2oo- Thus, our temperature estimates extracted within 
a 1 Mpc projected radius will slightly overestimate the 
full cluster temperature by ~ 10-20% relative to those 
derived in the sim ulations to select Coma-like clusters 
(|Bialek et al.l 12001). Therefore, the T^ values listed in 
column 6 of Table Q] need to be decreased by 10-20% 
to be accurately compared to the simulated Coma Clus- 
ter progenitor temperatures mentioned in the previous 
section. 
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3.1.2. X-ray Point Sources (XPSs) 

XPSs were identified using the C1AO tool wavdetect 
with a threshold limit set to 10 -6 which corresponds to 
one false source detection per 10 6 pixels (roughly the 
number of pixels in an unbinned image of one ACIS chip) . 
Each potential source above this threshold was scruti- 
nized individually. Broadband (0.3-8.0 keV) counts were 
extracted within a 95% encircled energy radius estimated 
for a monochromatic energy source of 1.5 keV at the ob- 
served off-axis angle. Background annuli regions were 
2.0 times the extraction radius. This background sub- 
traction does not account for the slow change of T x with 
radius within a cluster. However, a more sophisticated 
background subtraction using annuli at the XPS clus- 
ter radius show resulting differences less than the quoted 
photon statistics 

Conversion from net counts to flux was estimated for 
the XPSs using XSPEC assuming a power-law spec- 
trum with photon index of T=1.7 (Ng oc E~ r ). To 
estimate the X-ray flux limit for each observation (Col- 
umn 7 of Table [T]) , we estimated the point source flux 
for a 3cr detection above the background noise near 
the survey edge (R = 1 Mpc). Abell 2111 has a for- 
mal 3tr detection limit of £o.3-8.ofceV ~ 1-3 x 10 42 
ergs s _1 , slightly greater than our stated survey threshold 
of L .3-8.ofceV > 10 42 er g s s_1 - However, we feel com- 
fortable leaving this cluster in our sample because the 
CIAO wavdetect routine identified a point source with 
Lo.3-s.okeV ~ 8 x 10 41 ergs s _1 if the source is located 
at the clu ster redshift. P oissonian errors were calculated 
using the iGehrelsl (|1986l ) approximation. X-ray sources 
with SNR>3 and with expected K-corrected, rest-frame 
^0.3-8.0/ceV > 1 x 10 42 ergs s _1 are catalogued as po- 
tential cluster AGN candidates. We note that standard 
detection methods described above may miss XPSs asso- 
ciated with the BCG due to misidentifying AGN emis- 
sion as unresolved "cool core" ICM emission; i.e., the 
cusp of the ICM emission can mask the presence of an 
AGN. Therefore, we can make no solid claims about X- 
ray AGN in the bulk of the individual BCGs (see §g]for 
a further discussion.) 

3.2. Radio Imaging 

Several of our low-z clusters fields have been surveyed 
by the Faint Imag es of the Radio Sky at 20cm (FIRST; 
iBecker et al l [19951) and / or the NRAO VLA Sky Survey 
(NVSS: ICondon et "all 11998( 1 to detect radio galaxies 
with P\aghz > 3 x 10 23 Wm~ 2 . For our redshift range, 
the resolution of the FIRST maps may underestimate 
the total radio power of galaxies if any extended flux 
is close to the noise limit of the maps. Conversely, two 
radio sources with small angular separation may be cata- 
logued as one object by FIRST. To mitigate these issues, 
we obtained additional 1.4 GHz VLA maps (Program ID 
AS873) and re-analyzed 1.4 GHz VLA archival observa- 
tions obtained in its A-configuration. 

In September/December 2006 we obtained VLA con- 
tinuum observations at 1.4 GHz (50 MHz bandwidth) 
for four clusters (MS 0440.5+0204, Abell 2111, MS 
1455.0+2232, and Abell 1995). Abell 1995 was observed 
for 75 minutes in B-array, while the remaining three low- 
z clusters were observed for 30 minutes in C-array. The 
duration of a typical target scan was 15 minutes, brack- 



eted by 1 minute scans of a nearby phase calibrator. A 
flux calibrator (e.g. 3C286) was observed for 3 minutes 
at the beginning and end of each observing session. 

We used the NRAO Astronomical Imaging Processing 
System (AIPS), version 31DEC07, in the usual manner 
to flag, calibrate, transform and clean the images, so that 
radio sources could be detected manually. Source flux 
densities are estimated with the AIPS task TVWIN and 
IMEAN. The typical 3<r limit of our observations and 
FIRST is 0.3-0.5 mJy and 0.4 mJy, respectively. For 
Abell 963, the 3cr limit of the FIRST image is higher 
(~ 1.1 mJy) due to a bright background radio source 
with S\aghz ~ 1.4 Jy, but the observations still re- 
sult in a radio power limit below our required threshold. 
Higher resolution A-array maps for three clusters (MS 
0440.5+0204, MS 2137.3-2353, Abell 370) were reduced 
in a similar manner. In some instances, these higher reso- 
lution maps reveal two distinct radio sources for a single 
FIRST detection. While these higher resolution maps 
were used to provide the best optical identifications and 
to resolve close pairs of sources, the lower resolution maps 
(FIRST or other B/C-configuration maps) were used for 
flux determinations. 

In Table [T] columns 10-11 list the 1.4 GHz flux density 
and K-corrected luminosity limits for the individual clus- 
ters at the edge of our 1 Mpc survey region (K-correction 
assumes F v oc v~ a and a=0.7). These limits allow de- 
tection of many lower radio power FR 1 sources while 
excluding the lower luminosity radio sou rces due to star 
form ation (Li ±ghz < 10 22 - 75 W Hz" 1 : [Morrison et all 
[2001 . Given the very low radio power limits in col- 
umn 11, Table [TJ we do not expect that we have missed 
sources due to the presence of strong point sources in or 
near the primary beam, partially resolving the source or 
to beam-smearing in the outer regions of the fields, ex- 
cepting MS1358. 4 +6245. Due to a str ong, nearby radio 
source in this field, IStocke et al.l (|1999h only detected the 
BCG at 1.4 GHz. This 1.4 GHz map has a larger radio 
power limit (PiaGHz > 9.7 x 10 23 W Hz -1 ) for potential 
cluster radio galaxies than the o ther clusters in our sam- 
ple. However, in a 5 GHz map IStocke et al.1 (| 1999( 1 did 
not detect any additional sources within 5' of the BCG 
down to a 3a limit of 0.2 mJy, which corresponds to 1.4 
GHz flux limit of 0.5 mJy and Y lAGHz > 1.7 x 10 23 W 
Hz" 1 . Thus, we include MS1358. 4+6245 in our cluster 
sample with only one detected radio source in the BCG. 

3.3. Optical Datasets 
3.3.1. Multi-band Imaging 

Two-color images and photometry are publicly avail- 
able from the the Sloan Digital Sky Survey Data Re- 
lease 6 (SDSS DR6: IAdelman-McCarthv"et~all[20M) , the 
Canadia n Network for Ob servational Cosmology images 
(CNOC; lYee et al.1 fl996T l . an d/or the Chandra Multi- 
wavelength Project (ChaMP; [Green et al.ll2004H for 10 
of our 11 clusters in this paper. 

For 6 of our 11 clusters, SDSS DR6 provides five color 
imaging (ugriz) and photometry. The SDSS limiting 
magnitude in (g,r,i) is (22.2,22.2,21.3), adequately deep 
to detect M>M* cluster gala xies at z=0.4. SPSS pho- 
tome tric uncertainty is 1% ( Adel man-McCarthv et all 
I2008H . The CNOC cluster redshift survey examined 16 
X-ray luminous clusters with 0.2<z<0.55. For 3 of our 
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11 clusters, we utilized CNOC two-color photometric 
(Gunn g,r) and spectroscopic catalogs that are based 
on MOS observations from the Canada- Fra nce-Hawaii 
Teles cope (CFHT) between 1993-1994 (see lYee et al.l 
[19961 for details). For z<0.45 the CNOC limiting 
Gunn-r magnitude reaches down to M*+l with an 
individual source magnitude uncertainty of ~ 0.3 
mag. ChaMP is a wide-area survey of serendipitous 
X-ray sources in >100 archival Chandra observations. 
For MS 2137.3-2353, we use ChaMP's multi-band 
optical images (Sloan g,r,i), obtained at NOAO 4 m 
with the Mosaic CCD cameras and publicly available at 
|http://hea- www.cfa.harvard.edu/CHAMP /CHAMP.html| 

In November 2007 we obtained Sloan (g,r,i) images of 
Abell 370 with SPIcam on the Astrophysical Research 
Consortium (ARC) 3.5m Telescope at Apache Point Ob- 
servatory (APO). The FOV of SPIcam is 4.8 arcmin 2 ; 
therefore, we required multiple pointings of Abell 370 
to image the entire 1 Mpc radius region (FOV ~ 6.5 
arcmin 2 ). We used NOAO's Image Reduction and Anal- 
ysis Facility (IRAF), v2.14.1, in the usual manner to bias 
subtract and flat field individual images. Images were 
aligned and stacked into a large mosaic image with the 
IRAF package MSCRED. The limiting r-band magni- 
tude (~ 22.4) reaches down to M*+l with an individual 
source magnitude uncertainty of ~ 0.3 mag. For APO 
and ChaMP images, w e used the Picture Processing Pro- 
gram (PPP; lYedll99ll ) to detect and classify objects into 
stars and galaxies. PPP is a robust detection algorithm, 
especially in crowded fields, such as galaxy clusters. Ad- 
ditionally, total source magnitudes are calculated using 
optimal extraction radii that are based on photometric 
growth curves. 

3.3.2. Spectroscopy 

SDSS and/or CNOC spectroscopic data are available 
for several candidate cluster XPSs and radio galaxies. 
The SDSS spectral coverage is 3800-9200 A with a res- 
olution of 1800-2200. The SDSS spectroscopic pipeline 
(spectrolD) cross-correlates individual wavelength and 
flux-calibrated spectra with various template spectra 
(stellar, quasar, and emission-line galaxies) to determine 
source re dshifts via emission line s and/or absorption fea- 
tures fsee lStoughton et al. I [20021 for details). The CNOC 
spectral coverage is 4650-6100 A with a dispersion of 
~ 3.45A/pixel, yielding a spectral resolution of 16.5 A 
with a 1.5" slit width. The CNOC survey also incor- 
porates cross-correlation techniques to determine source 
redshifts that ha ve velocity unce rtainties between 100- 
130 km s" 1 fsee lYee et all 119961 for details). In addi- 
tion to these two surveys, we searched the literature (via 
NASA's Astrophysics Data System Bibliographic Ser- 
vices) and extragalactic databases (e.g. NASA/IPAC 
Extragalactic Database) for published redshifts and/or 
spectra of our potential cluster targets. 

For remaining cluster AGN candidates without pub- 
lished redshift and/or spectra, spectroscopic observations 
were obtained for all objects with Sloan r < 20.8 using 
the Double Imaging Spectrograph (DIS-II) at the APO 
3.5m between 2006-2009. Observations were obtained in 
its low resolution (~ 6 A) mode which provides complete 
spectral coverage from 3700-9000 A. Comparison spec- 
tra (He, Ne, Ar) were routinely obtained between targets 
and/or long exposure sets and subsequently were used to 



calibrate the wavelength scale of science spectra. Expo- 
sure times ranged from 45 minutes to 2.5 hours, achiev- 
ing a SNR > 5 per pixel in the science spectra. Spec- 
troscopic images were reduced with IRAF in the usual 
manner (bias subtraction and flat fielding) and source 
spectra were extracted and calibrated with the IRAF 
APEXTRACT package. 

Source redshifts were calculated from individual wave- 
length measurements of typical absorption features (e.g. 
the 4000A Ca break, G-band, Mg lb, Na I) and/or 
emission-line features (e.g. [O II], [O III], Ha). The 
radial velocity accuracy from a typical spectral line is 
~ 300 km s _1 (or Az ~ 0.001). At this writing, spectra 
for radio galaxies and XPSs are complete to M r < —20.8 
(M * 20.8 for typic al SDSS DR4 rich cluster galax- 
ies; [Hansen^FaDIiOOil). We define passive galaxy spec- 
tra as those showing only stellar absorption lines (e.g., 
the 4000A Ca break, Mg lb, Na I), while active galaxy 
spectra have detectable emission lines (e.g., [O II], Ha, 
[O III]). We do not find any example of a starburst op- 
tical spectrum among either the radio galaxies or XPSs 
as determined by comparisons of strong emission-line lu- 
minosities (e.g. [0 III] 5007 A comparable to H/3; see 
iStocke et al.lll99lt ). This method is indicative, and not as 
robust as "BPT diagnostic plots" ([Baldwin et al.lll981h 
that utilize several emission-line features to separate star- 
burst galaxies from AGN. 

3.3.3. Cluster Red Sequence Galaxies 

The color-magnitude diagrams (CMDs) of clusters gen- 
erally display a tight color sequence for the cluster el- 
lipticals, referred to as the cluster red sequence (CRS). 
We estimate the number of CRS galaxies within 1 Mpc 
of the cluster center and with M r < —20.8 (where our 
spectroscopy becomes incomplete) in the following man- 
ner. Foreground and background galaxies are statisti- 
cally subtracted from e ach cluster, u s ing S loan r-band 
background counts from lYasuda et all (|2001h . We deter- 
mine the mean CRS color by fitting a bimodal Gaussian 
distribution to both the red and blue galaxy populations 
as a function of color. Using only galaxies with colors 
within 3cr of the mean CRS color, we define the CRS 
more precisely by fitting a line to the observed color and 
magnitude of these sources and thus verifying the well- 
known " tilt " of the CRS for cluster color-m agnitude dia- 
grams (e.g- JVisvanathan fc Sandagelll977[ ). For galaxies 
with M r < —20.8, the large majority of blue galaxy con- 
taminators are foreground galaxies and not the fainter 
blue cloud galaxies with M r > M* . The estimated num- 
ber of CRS galaxies will be used to compute AGN frac- 
tions described in § 14.11 

4. RESULTS 

Within 1 Mpc of the cluster X-ray emission centroid, 
we find 8 XPSs with £0.3-8. okeV > 1 x 10 42 ergs s" 1 , 
confirmed to be cluster members within the 3<r velocity 
dispersion of the cluster (from the literature) and M r < 
-20.8, with the exception of Abell 963 X2 with M r = 
—20.3, which we include despite being 0.5 mag fainter. 
We also find 20 radio galaxies with P\aghz > 3 x 10 23 
W Hz -1 and M r < —20.8, of which six are located in the 
BCG. 

Tables [2H3] lists the basic data for the XPSs and ra- 
dio galaxies, respectively, in our sample. The sources 
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are identified by cluster name and a consecutive number- 
ing which includes foreground and background sources. 
In Table the letter identifiers in parentheses are 
used in subsequent plots and discussions for easy cross- 
referencing. Both tables include the spectroscopic red- 
shift (mostly from this work), observed magnitude of the 
host galaxies in which they are found, the projected ra- 
dial distance from the cluster center as defined by the 
diffuse ICM X-ray emission peak in h^ 1 Mpc (called " ra- 
dius"), and observed broadband (0.3-8 keV) X-ray lumi- 
nosity or limit. Table [2] also includes the net broadband 
(0.3-8.0 keV) counts for each XPS. 

Table [3] also includes the radio source flux density at 
1.4 GHz and K-corrected 1.4 GHz radio luminosity. Host 
galaxy colors and cluster "red sequence" (CRS) colors 
(see § I3.3.3j) were obtained using Sloan (g-r) for z<0.3 
clusters or (r-i) for z=0.3-0.4 clusters (photometric cal- 
ibration errors are typically 2% in these colors). For the 
3 EMSS clusters, CNOC photometry provides Gunn (g- 
r) colors. We note that MS 0440.5+0204 R3 & R4 and 
Abell 370 Rl & R3 are examples of two distinct sources 
detected in the higher resolution VLA maps, but were 
identified as a single source in FIRST images. 

For radio galaxies with detectable X-ray emission, the 
X-ray luminosities (Table[3j column 12) were determined 
in the same manner as for the detected XPSs described in 
§ 13.1.21 The upper X-ray luminosity limits for the BCGs 
and radio galaxies with X-ray non-detections above our 
limits are estimates as 3er above the ICM emission at 
their locations. For Abell 1758, the ICM X-ray emis- 
sion is both diffuse, allowing a good estimate for the 
upper limit on the X-ray luminosity of the BCG, and 
double-peaked, indicative of a merging system. In this 
case we have determined the projected distance of the 
radio galaxies to be from the nearest diffuse ICM emis- 
sion center. Thus, A1758-R4 is identified with the BCG 
in the NW X-ray clump and A1758-R6 is 2/3 Mpc from 
the other clump center to the SE. Otherwise, there are 
no ambiguous cases in the projected distance (in Mpc) 
from the center of the ICM X-ray emission. The radio 
upper limits for the XPS locations are set at 3er above 
the background noise and so are conservative. 

Due to the difficulty of detecting XPS at the center 
of the diffuse X-ray emission, we probably have under- 
estimated the number of cluster XPSs. We did detect 
one BCG in both radio and X-rays, Abell 1758 R4/X1 - 
listed in Tables [2H31 where the ICM emission is particu- 
larly diffuse. Since 6 of 11 BCGs were detected as radio 
galaxies (refer to Table [3]), we might expect that a simi- 
lar fraction of BCGs are XPSs. However, since our radio 
and X-ray detection limits allow the detection of only 
half the total number of XPSs as radio galaxies (refer to 
Table [2]) we conservatively estimate that ~ 3 BCGs are 
XPSs at £o.3-8.ofceV > 10 42 ergs s _1 in our sample, one 
of which we have already identified (A1758 R4/X1). 

Figure [2] displays a composite color-magnitude dia- 
gram for our eleven low-z clusters. All but one ra- 
dio galaxy lies on the CRS. The exception is the BCG 
in MS 1455.0+2232 for which extremely luminous, ex- 
tended [O II] 3727 A emission is present in the g-band, 
making it appear to be somewhat bluer than the CRS 
(jDonahue et a.1.1 fl992f ) . Also excepting one case (Abell 
963 XI, object A in Figure [2]), XPSs are hosted by early- 
type galaxies on or near the CRS. Abell 963 XI has 



an optical spectrum consistent with a Seyfert nucleus 
(with emission-line luminosity of [OIII]>>H / 3). Given 
the presence of only one Seyfert projected within 1 Mpc 
of a cluster core and ~ 745 km s _1 relative radial veloc- 
ity for this Seyfert ( Abell 963 cluster veloci ty dispersion 
of cr„ = 1350 km s" 1 ; iStruble fc RflodHl99flh . this source 
is consistent with being a cluster member either fore- 
ground or background to the core at a physical radius 
> 1 Mpc. A second XPS (Abell 963 X2, object D in 
Figure [5]), the least luminous XPS host galaxy in our 
sample, is several tenths of a magnitude bluer than the 
red sequence, but, nonetheless, possesses a passive ab- 
sorption line spectrum. This galaxy will be discussed in 
detail in § [431 

4.1. AGN Fraction and Radial Distribution of Cluster 
Red Sequence (CRS) Radio Galaxies and X-ray 
Point Sources 

Table [4] lists our estimates for the number of 
background-subtracted CRS galaxies for each cluster (see 
^ 13.3.31 for details), the total number of cluster radio and 
X-ray sources, the radio and X-ray source fraction, color 
used to determine the CRS number, and the optical sur- 
vey utilized. Background-subtracted CRS galaxies num- 
ber ~ 665 in our entire cluster sample. The average AGN 
fraction for CRS galaxies in our sample is 3% for radio- 
selected sources and 1% for X-ray selected sources within 
1 Mpc of the cluster center. 

The cumulative radial distributions of the cluster radio 
galaxies and XPSs, as shown in Figure [3j reveal these 
sources to be more centrally-concentrated in the clus- 
ter relative to CRS galaxies as a whole. Eighty percent 
of radio galaxies and 60% of XPSs are located within 
500 kpc of the cluster centers and appear more cen- 
trally concentrated than the CRS galaxies even in that 
region. A luminosity weighting of the CRS distribution 
using the bivariate rad io luminosity function derived by 
lAuriemma et al.l (|1977h is indistinguishable from the un- 
weighted CRS galaxy distribution; i.e., there is no evi- 
dence for significant mass segregation in these clusters. 
Thus, the significant difference between the radio galaxy 
+ XPS population projected radial distances and the 
weighted CRS galaxy radial distances remains as shown 
in Figure [3] 

A two-sided Kolmogorov-Smirnov (K-S) test between 
the full X-ray + radio population and the CRS galaxy 
population as a whole are inconsistent with being drawn 
from the same parent population at > 99.999% confi- 
dence level (K-S D-statistic = 0.36 and probability = 
3.3 xl0~ 6 ). The K-S probability is even lower if we only 
examine the cluster radio sources (K-S D-statistic = 0.44 
and probability = 4.5 xlO -9 ). The K-S probability is no- 
ticeably larger for just the cluster XPSs (K-S D-statistic 
= 0.26 and probability = 1.9xl0~ 3 ), but this compari- 
son is still significant at the 99.8% confidence level. To 
account for the likely presence of XPSs in BCGs which 
are not detected due to the presence of an ICM "cool 
core" cusp in the diffuse X-ray emission, we add two 
XPSs at r=0. The K-S probability for this augmented 
XPS distribution (K-S D-statistic = 0.38 and probabil- 
ity = 7.1 xl0~ 7 ) is comparable to the original X-ray + 
radio population and only strengthens our conclusion of 
a centrally concentrated population of XPSs. 

Within 500 kpc 16 radio galaxies in our sample make 
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up 6±2% of the bright (L> L*) ~ 265 CRS galaxies 
(40% of 665), while the 4 XPSs (excluding Abell 963 
XI) make up 2±1%. The combined radio and XPS pop- 
ulation make up 8±2% of the bright (L> L*) ~ 265 CRS 
galaxies. Our result strongly suggests that the triggering 
of a radio-loud AGN or an XPS in these cluster galaxies 
is due to some, as yet undetermined, mechanism related 
to the ICM density s uch as in the Bondi accretion model 
of lAllen et, all (f2006T ). 

4.2. X-ray Properties 

Of the 8 detected XPSs (Table [2]), we detect 7 XPSs 
with 25 — 115 net broadband (0.3-8.0 keV) counts and 
one XPS (Abell 963 XI) with - 2150 net counts. This 
XPS is the bright (M r ~-22.8), blue (almost 1 magni- 
tude bluer than the CRS) galaxy in Figure [2] (Object A) 
and possesses a typical Seyfert optical spectrum. Using 
XSPEC, the spectrum of Abell 963 XI is well-modeled by 
a power-law with spectral index T = 2.0 ± 0.1, emission 
from the 6.4 keV Fe K-a line with an equivalent width of 
1.7 keV, and no intrinsic absorption. But this spectrum 
also can be fit by a Raymond-Smith thermal model with 
kT = 3A±° Q ;| keV, but at a higher reduced-x 2 value (1.18 
vs. 1.03). Given the clear Seyfert-like optical spectrum 
and likely power-law X-ray spectrum we identify Abell 
963 XI as a clear AGN. 

But due to low photon counts in the other 7 XPSs, 
we instead compute X-ray colors to estimate their spec- 
tral slope and amount of intrinsic absorption. Figure [4] 
displays the ratio of two X-ray colors. The three bands 
(SI, S2, and H ) are defined for similar energy ranges 
as in lKim et all (|2004h and lMartini et alj (|2006f) . Over- 
laid are the expected colors for a power-law spectrum 
with spectral slopes of T = 1.0 — 4.0 and intrinsic 
riH = 0, 10 21 , 10 22 , 10 23 cm -2 . Each individual source 
is displayed and labeled as in Figure [2] and in Table [2 
the composite X-ray color (filled square) is also displayed 
using all the XPSs without detectable radio emission and 
the XPSs with detectable radio emission (filled triangle), 
excepting Abell 963 XI (Object A). With the exception 
of MS0440 X8 (Object G) and RXJ0952 X6 (Object H), 
the XPSs are cons istent with a typic al AGN photon in- 
dex of T - 1.75 (jTozzi et al.|[200l ). MS0440 X8 has 
relatively few hard (2.5-8.0 keV) counts and has consis- 
tent colors for source with r ~ 3, but as with the others, 
little evidence for obscuration in the X-rays. RXJ0952 
X6 has approximately equal counts in the soft (0.5-2.5 
keV) and hard (2.5-8.0 keV) bands and has consistent 
colors for a source with r ~ 1. 

In addition to the analysis above, a composite X-ray 
spectrum in five discrete rest-frame bands is consistent 
with a power-law model with T = 1.6 ± 0.2 with no ev- 
idence for intrinsic absorption. Thus, if these XPSs are 
AGN, they are unob s cured AGN . There is no evidence for 
X-ray spectral differences for XPSs with and without ra- 
dio detections, and, although the statistics are modest, 
this evidence also weighs in favor of all of these XPSs 
being AGN of some sort. 

4.3. X-ray to Radio Luminosity Ratios 

We find that XPSs and radio galaxies rarely overlap 
in our sample. There are only 3 XPSs which are radio- 
detected, suggesting two populations. Are we really see- 
ing two populations of sources or is this just a function 



of our flux limits? To investigate this further, we deter- 
mined the X-ray to Radio Luminosity ratios for our clus- 
ter sources (Figure [5]) an d compared them to op t ical e l- 
liptical "core" galaxies in IBalmaverde fe Capetti| {2 0061. 
more luminous 3C/FR 1 radio galaxies in iDonato et al.l 
(l2004fl and Seyferts hosted by early-type galaxies in 
iCapetti fe Balmaverdei (|2007[ ). A loose linear relation- 
ship exists in Figure [5] for the low (optical "core" 
galaxies) and higher p ower (FR Is) ellipticals of the 
IBalmaverde fe Capettif (|2006l ) and IDonato et all 1)20041 ) 
samples, respectively. The X-ray limits on the radio 
galaxies in our sample are consistent with other FR Is; 
i.e., we would require perhaps a factor of ~ 3 improve- 
ment in X-ray sensitivity before we would expect to de- 
tect all of these radio galaxies. 

On the other hand, the XPSs in our sample would be 
easily detected in the radio if they were similar to FR 1 
radio galaxies; indeed, they should have been detected at 
power levels comparable to or significantly higher than 
the radio galaxy detections we did make. So, we can rule 
out that the radio galaxies and the XPSs are drawn from 
the same parent population. But the XPSs do have X- 
ray luminosities and radio luminosity limits that place 
them among the Seyferts in Figure [5] If these XPSs 
are some sort of Seyfert galaxy, an order of magnitude 
improvement in our radio flux limit probably would be 
required to detect all of them. Indeed, if the XPSs are all 
Seyferts, it is somewhat surprising that we did not detect 
some of them already given the distribution of Seyferts 
in Figure El 

Are the XPSs simply radio-quiet AGN, similar to 
Seyfert galaxies? While they possess X-ray luminosities 
of Seyferts, they do not possess the emission-line spectra 
of Seyferts nor is there unambiguous evidence for ob- 
scuration. Even if we use the X-ray to Ha flux ratios 
(mean value = 7.3) of much lowe r luminosity AGN like 
in the nearby sample of lHol (|2008l ). we would expect Ha 
and other emission line fluxes for our sample of XPSs at 
L_Ha > 10 41 ergs s -1 , which we easily exclude using our 
3.5m spectra. We have also gone through the exercise of 
determining if the Seyfert spectrum (i.e., blue power-law 
and luminous emission lines of [OIII] and Ha) of Abell 
963 XI would be detectable if scaled down by the dif- 
ference in X-ray luminosity between this sources and the 
others (using standard scaling relations observed for X- 
ray emitting Seyferts). We find that a Seyfert-like AGN 
of power 10% of Abell 963 XI would be detected eas- 
ily in our optical spectra and a 3% power-level probably 
would be detected as well. But it might be possible to 
"hide" a Seyfert 2 type spectrum (line luminosities can 
be an order of magnitude lower than many Seyfert Is), 
which leads to the suggestion that these could be highly- 
obscured AGN. However, as we showed above, there is 
no strong evidence for obscuration in these XPSs based 
on their X-ray spectra. Nor are their optical broad-band 
colors significantly redder than the CRS although a low 
luminosity AGN with M r ~ — 18 could be hidden in these 
luminous galaxies without affecting the broad-band col- 
ors within our photometric errors. We conclude that the 
XPSs are unlikely to be similar to Seyferts or that they 
are highly-obscured AGN of any description. 

4.4. Comparison to Other Cluster AGN Studies 
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How do our results compare to other cluster AGN 
studies? Cluster radio g alaxies have been kn own to re- 
side in passive galaxies (|Morrison et al.l 12003). thus our 
findings on their host galaxies are not surprising. Sta- 
tistical studies of cluster radio galaxies have also re- 
vealed enhanced concentration towards cluster centers 
(e.g., iLedlow fc Owenl fl995t iLin fc Mohrl I2001 similar 
t o what we have found. 

ILedlow fc O^enl ((19951 ) surveyed ~300 low-z (z<0.09) 
Abell clusters to find radio galaxies within ~0.6 Mpc 
of the cluster optical center. The authors find that 
the surface density distribution of radio galaxies with 
log(Pi.4G#.z)>23.23 is more centrally concentrated than 
an expected King mod el distribution of cluster galaxies. 
ILedlow fc Owenl Tl995) argue that the central concentra- 
tion of radio galaxies in clusters is a direct result of (1) 
clusters having luminous host galaxies residing in their 
cores and (2) more luminous galaxies having a higher 
proba bility of hosting a radio source (jLedlow fc Owenl 
1996). Our radio galaxy AGN fraction within 0.5 Mpc 
of the cluster center is <~6%, consistent with their esti- 
mated fraction (~8%) o f elliptical galaxies hos ting radio 
s ources (see Figure 6 in lLedlow fc Owenlll995l ). 

ILin fc Mohrl (|2007h surveyed 573 X-ray detected clus- 
ters that were also observed with NVSS. The authors 
also find the surface density distribution of radio galax- 
ies (Piaghz)> 10 23 W Hz -1 ) within 5r 20 o to be more 
centrally con centrated than Tw o Micron All-Sky Sur- 
vey (2MASS; Uarrett et aTl l2000f ) K-band galaxies with 
Mk<-24. They estimate the radio-active fraction (RAF) 
of cluster and field galaxies with Piaghz > 10 23 W Hz -1 
and M K <-24 to be 4.9 ±0.5% and 1.3±0.5 %, respectively. 

Note for M* 20.8 (|Hansen et al.ll2005D an d (r-K)~2.9 

for SDSS early- type galaxies ((Chang et al.ll2006T ). then 



M^- ~ —23.7. These authors surveyed a larger region 
around clusters (R~ r2oo) that could account for the 
smaller cluster galaxy RAF compared to our value. Also, 
there is both a luminosity and a Hubble type bias to this 
difference since many field galaxies are late-type systems 
which do not harbor radio galaxies. 

Studies of X-ray point sources in the fields of small 
samples of galaxy clusters (e.g., ICappi et al.l 120011 : 
iBranchesi et al.l I2007D revealed a slight excess ~ 2a in 
the exp ected surface density c ompa red to cluster-free 
regions. iRuderman fc Ebelingj ((20051 ) surveyed 51 clus- 
ters with 0.3<z<0.7 and L .5-2keV >10 43 ergs s -1 and 
found a large excess in the composite surface density 
of XPSs within 3.5 Mpc of the cluster center com- 
pared to their 20 control fields. In relaxed clusters, 
i.e. clu sters with smooth, symmet ric X-ray emission 
profiles, IRuderman fc Ebelingj ((2005T ) find a prominent 
excess within 0.5 Mpc, whereas clusters with disturbed 
X-ray morphologies appear to have their XPS excess dis- 
tributed mo re uniformly wi t hin th e larger 3.5 Mpc sur- 
vey region. iGilmour et al.l (|2009l ) analyzed 148 galaxy 
clusters with 0.1<z<0.9. Within a 1 Mpc projected re- 
gion they find a 3a excess in XPS counts compared to 44 
control fields, resulting in ~ 1 . 5 XP Ss per cluster across 
their sample. IGilmour et alj ([2009) find more AGN in 
the outer 0.5-1.0 Mpc regions than within 0.5 Mpc. In 
these larger statistical surveys all XPSs (after a statis- 
tical background subtraction) are assumed to be at the 
cluster redshift, but little or no redshift information ex- 
ists to verify cluster membership. 



However. (Martini et al.l (|2006l I2007D has obtained red- 
shifts for cluster AGN candidates. They examined eight 
low-z clusters (0.05<z<0.31) with cluster ICM tempera- 
tures ranging from 4 to 16 keV (including MS1008 which 
overlaps with our sample) and with 3 fields surveyed out 
to a 1 Mpc projected radius. Although their X-ray lu- 
minosity limit is one order of magnitude lower and their 
magnitude completeness limit is one magnitude fainter 
than MJj, their R-band images show that most of their 
X-ra y point sources are hosted by galaxies with colors 
(see (Martini et alj (20071 Fig. 3 ) consistent with (B - 
R) = 1 .57 ±0.2 of local e l lipticals (|Fj jkirgita_elalJ 119951) . 
So while iMartini et all (|2006l . [2007D do find a few blue 
cluster AGN with Seyfert-like spectra (2 out of 6 X-ray 
sources with Lx > 10 and located within 1 Mpc of 
the cluster cores), their numbers are small. If we assume 
that the XPSs in our sample are all AG N, then our X- 
ray AG N fraction of ~2% is consistent with lMartini et aTl 
(2006) for L0.3-s.0keV > 10 42 ergs s 1 . It is interesting 
to note that although we have a more robust selection 
method for our low-z cluster sample, t he radio galaxy 
AGN frac tion of [Cedlow fc Ow ej (|l995l) and XPS AGN 
fraction of lMartini et al.l (|2007| ) are similar to our values. 

5. THE NATURE OF CLUSTER X-RAY POINT 
SOURCES (XPSS) 

As stated in the previous section, our cluster XPSs 
do not appear to have optical nor X-ray signatures typ- 
ical of Seyferts. However, the XPSs have many of the 
properties of low luminosity BL Lac Objects of the 
class now called "High-e nergy-Peak" BL Lacs (HBL; 
iPadovani fc Giommilll995l; first called X-ray selected BL 
Lac Objects: IStocke et al.l (T985). With X-ray luminosi- 
ties ~ 10 42 ergs s _1 an HBL would be expected to have 
a radio luminosity of 10 21-23 W H z -1 based upon HBL s 
at higher total luminosities (e.g.. [Perlman et alj fl996). 
at or below our radio detection limits for this sample. 
Indeed, 3 of these XPSs are radio-detected. Most HBLs 
previously studied have X-ray and radio luminosities 2-3 
orders of magnitude more luminous than our detections 
and still only sometimes exhibit the featu reless power- 
law o ptical spectrum of classical BL Lacs ((Rector et al.1 
1999). Mostly, HBLs have optical spectra of passive ellip- 
tical galaxies with emission line luminosities < 10 40 ergs 
s " 1 , consistent with w hat we find for the cluster XPSs. 

iRector et al.1 ((19991 ) investigated the general proper- 
ties of low-luminosity HBLs, finding that their properties 
merge with normal ellipticals, with normal (i.e., ~ 50%) 
Ca I I H fc K "break strength s" seen in non-AGN ellipti- 
cals. lUrrv fc Padov ani ( 1995) reviewed the other proper- 
ties of these objects whi ch find them in giant (M r = -22 to 
-25) elliptical galaxies (|Wurtz et alj Il997t iFalomo et all 
1999); and in rich groups to mod erately rich cluster s 
at similar redshifts to our sample ((Wurtz et al.l I1997T ). 
HBLs show no evidence for internal X-ray obscuration 
and some show quite soft X-ray spectra like MS0440 X8 
(|Perlman et al.lll996l ). 

The detailed optical properties of Abell 963 X2 (source 
D in Figures [U & O also support the HBL classifica- 
tion of the XPSs. Figure [2] shows that this XPS host 
galaxy is ~ 0.5 mag bluer in (g-r) than its associated 
CRS. Further to the blue, the discrepancy enlarges with 
Sloan (u-g)=0.64±0.26 while the CRS galaxies in Abell 
963 have Sloan (u-g)=l. 4-2.4. Deeper (u,g,r) imaging 
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of Abell 963 X2, obtained by us at the APO 3.5m tele- 
scope, confirms its anomalous colors blueward of the Ca 
II break. The radio/XPS source MS1455 X2/R1 (optical 
object E in Figs. [2j [4j &[5| has colors slightly bluer than 
the CRS, and so could also possess non-thermal emis- 
sion in the blue. Table [5] compares the observed colors 
of these two "blue" XPSs to the expected colors for a 
cluster galaxy of similar r-band luminosity and redshift. 

This excess blue flux also is confirmed by our optical 
DIS spectrum which sh ows a slightly diminished C all 
HfcK "break strength" (|Dressler fc Shectmanl Il987f) of 
0.4 s imilar to other low-luminosity HBLs ([Rector et aTl 
fl999T) . Since this is the lowest luminosity host galaxy 
among our XPS population, this is the most likely galaxy 
for which the blue power-law continuum of a BL Lac 
would be expected to be visible. For the other XPSs in 
our sample, the Call H&K "break strength" ranges be- 
tween 0.4-0.5, in support of the low-luminosity BL Lac 
hypothesis. Confirmation of this hypothesis can be ob- 
tained by detecting optical polarization in the B-band 
(blueward of Ca II H&K at this redshift in Abell 963 
X2). Thus, the observed properties of the XPSs match 
those of HBLs the best of any known AGN class. 

The X-ray l umino sity f unctions obta i ned fo r HBLs by 
iMorris et at] (|1991h and IRector et all (|2000h show no 
signs of flattening down to L^^IO 43 5 ergs s _1 (although 
theoretical luminosity functions of bea med sources all 
have a flattening at low power levels; lUrrv fc S hafcr 
1984) and so numerous lower power HBLs would be ex- 
pected in bright cluster ellipticals. In Figure El except- 
ing the one Seyfert (Object A), the other XPSs are 1.5-2 
orders of magnitude more luminous in X-rays than the 
sequence of FR Is and compact core galaxies. Thus, 
the approximate "excess L x " of the XPSs in our sample 
above the FR 1 sequence suggests a Doppler-boosting of 
X-rays by factors of 3 0-100. Taking the simpli fied model 
of HBLs described in lUrrv fc PadovanH (jl995), our line- 
of-sight is located well off the radio beaming axis (thus 
no beamed radio emission is observed) but within the 
7 > 3 X-ray beam. In this case the maximum kinematic 
Doppler factor expected for the X-rays is S ~ 2j and the 
ratio of observed-to-intrinsic X-ray luminosity is S p where 
p=3+a for discrete emitting "blobs" and p=2+a for a 
continuously emitting structure (where a is the intrin- 
sic spectral index in energy units for the source; here we 
conservatively assume a= 1). In the simplifying assump- 
tions of the HBL model we obtain a maximum luminosity 
boosting of 200-1300 depending upon the physical struc- 
ture of the emitting region (or even higher if T is greater 
than as estimated from the inferred open ing angles of 
25-60 degrees based upon source counts; IRector et al.l 
120001) . Since the most luminous HBLs have L x ~ 10 44 " 45 
ergs s _1 , the unbeamed X-ray emission is estimated to 
be 10 40-41 ergs s _1 , or even lower depending upon how 
these quantities scale with source luminosity. Thus, our 
estimate above that the XPSs are boosted in their X-ray 
emission by 30-100 times seems plausible. 

The most luminous HBLs are found in large-solid- 
angle "serendipitous" X-ray surveys and so are quite rare 
(~ 10~ 7 sources Mpc~ 3 at L x > 10 44 ergs s~ 4 from 
the IRector etail ((2000) XLF). A power-law extrapola- 
tion of that XLF yields ~ 10~ 5 BL Lacs per Mpc -3 
at L x > 10 42 ergs s~ x . Since the elliptical galaxy LF 
finds - 8 x 10~ 4 Mpc- 3 for L>L* (jMarzke et al.lll99l 



Urn et al.lll996l) , the expectation from observed space den- 
sities is ~ 1 BL Lac in 80 > L* elliptical galaxies. In 
this survey we have observed ~ 665 bright cluster el- 
lipticals and found 7 XPSs (not include Abell 963 XI) 
which translates to ~ 1 XPS per 95 bright ellipticals, re- 
markably identical to the expectations given the rather 
large uncertaint ies due to counting statistics in both our 
sample and the IRector et al.l (|2000f ) HBL sample. 

In conclusion, these XPSs possess all of the qualities 
(including space densities) expected for low luminosity 
HBLs and thus are related to the higher power radio 
galaxies in our sample in that both classes are radio- 
loud and possess relativistic jets that can carry energy 
and material into the surrounding ICM. If the identifi- 
cation of these XPSs as HBLs is correct, then both the 
radio-selected and the X-ray selected AGN in clusters are 
radio-loud sources, and there is no need to create a new 
category of AGN for similar sources seen in other sur- 
veys (e. g., deep XMM and Chandra serendipitous sourc e 
survevsT lBarger. Cowie. Mushotzky. fc Ric hards 2001b). 

A final known category of X-ray emitter into which the 
XPSs might be placed is an a ccreted "cool cor e" thermal 
corona in galaxy-sized halos (|Sun et al.l l2007). Thermal 
X-ray coronae of early type galaxies hav e been known 
since the X-ray observations of Einstein (|Forman et al.l 
1985). Examples of t his phenomenon have been seen 
in ne arby clusters (e.g.. lYamasaki et al.ll2002HSun et al.l 
2005), and not necessarily in the BCG. We expected to 
spatially resolve this type of emitting region since the 
nearest XPSs in our sample are at a distances where 
the 1 arcsec Chandra resolution is ~ 3 kpc. But, the 
surrounding diffuse hot ICM emission may make these 
cores more difficult to resolve at higher redshift, espe- 
cially when they are detected at such low count levels, 
so this test is not definitive for our data. Additionally, 
simulations of cool cores in in-falling cluster ellipticals 
are 1-2 orders of magn itude less than the X-ray emis- 
sion in these XPSs (e.g.. lVikhlinin et al.ll200H ISun et all 
2007). And while the X-ray spectra of the XPSs we have 
found can be fit by thermal models with T^ less than 
the surrounding medium, a power-law fit is better in all 
individual cases and composites (see § 14. 2|) . The typical 
observational signature for soft thermal coronae is the 
presence of the iron L-shell b ump at 0-8-1 keV , resulting 
m a poor power law fit fe.g- lSun et all[2007h . The one 
possibility for which we can imagine creating cool cores 
as luminous as >10 42 ergs s _1 is a remnant core of a 
group BCG that is in the process of falling into the clus- 
ter. But none of the XPSs we detect (except one cluster 
BCG which also has radio emission) is in a very bright 
CRS galaxy. We conclude that it is unlikely that these 
sources are remnant cool cores and that the best model 
for XPSs is HBL-type BL Lac Objects. 

6. DISCUSSION 

Based upon the observed properties of the 20 radio 
galaxies and the 7 XPSs (total of 24 AGN counting 
overlaps and excluding Abell 963 XI) in these 11 low- 
z clusters, we conclude that all these sources are radio- 
loud AGN. The radio galaxies have properties consistent 
with high luminosity FR 1 sources while the XPSs with- 
out radio emission are most consistent with being low 
power FR Is with X-ray emission that has been Doppler- 
boosted above our detection threshold. Thus, there is ob- 
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servational evidence that both AGN classes (X-ray and 
radio selected) possess relativistic jets that can heat the 
surrounding ICM through jet interaction. 

Since both classes have luminosity function determina- 
tions from previous work, we can extrapolate downward 
in both X-ray and radio flux below our threshold using 
previous data. We have normalized these extrapolations 
using our own sample of ~ 265 L > L* cluster ellipticals 
at projected radial distances < 500 kpc from the cluster 
ICM X-ray emission centroid, of which 6% (16 objects) 
are radio galaxies at I^xaGHz >3 X 10 23 W Hz™ 1 and 1% 
(3, excluding A963 XI) are XPSs at L x > 10 42 e r gs s" 1 . 

Based on Figure 4 of iMachalski fc GodlowskH (gOOO), 
we integrate the AGN radio luminosity function (RLF) 
for log(P)>23.4 (high power) and for 21.4<log(P)<23.4 
(low power) and expect 13 times as many low power ra- 
dio galaxies as high power radio sources in our sample 
(or ^210 low power radio sources). Based upon these 
values, we predict that ~ 85% of all bright CRS galax- 
ies within 500 kpc of cluster center possess radio sources 
with L 1AGHz > 10 21 - 4 W Hz" 1 . The XLF of HBLs pre- 
dicts even a higher number, all (100%) bright cluster 
ellipticals should possess an XPS at L T > 10 40 ergs s" 1 



But t he XLF of (Morris et all (|1991[ ) and iRector et all 
(2000) extends at constant logarithmic slope down only 
to 10 ergs s _1 . Since flattening of this function is ex - 
pected at or below 10 43 ergs s" 1 (jUrrv fc Shaferlll984f ). 
the number of predicted HBLs would be somewhat less 
than 100% and so in good agreement with the extrapo- 
lation of the RLF. Thus, the two extrapolations roughly 
agree in predicting that almost all bright CRS galaxies 
(> L*) at projected radial distance < 500 kpc contain 
radio-loud and X-ray loud AGN at h 1AGHz > 10 214 W 
Hz" 1 and L x > 10 40 ergs s _1 , respectively (see Fig. [5]). 
This agreement is additional evidence that our identifi- 
cation of the XPSs as low luminosity BL Lacs is correct. 
This result suggests that deeper X-ray surveys of nearby 
clusters will discover large numbers of XPSs in bright el- 
lipticals and very deep Chandra i maging of the Perse us 
Cluster appears to show just that (jSantra et al.l l2007). 

What impact do all of these AGN have on cluster heat- 
ing models? They certainly can provide a distributed 
source of heat which moves about the inner 0.5 Mpc or 
more of these clusters on a characteristic crossing time 
of a few x 10 8 yrs. But is their input miniscule com- 
pared with the AGN in the BCG? To evaluate the heating 
by these AGN we have assumed results from two recent 
studies of th e relationship b etween r adio luminosi t y and 
jet p ower bv lBicknelll (|1995f ) and bv lBirzan et~aT1 (|2004l . 
2008). First we assume that the magnetic field strength 
in the inner regions of these clusters is relatively constant 
so that the ICM heating rate is proportional only to the 
AGN radio luminosit y (i.e., a const ant efficiency factor 
for heating the ICM; lBicknelll [l995). This assumption 
allows us to use the observed RLF to predict that all 
cluster AGN below Li AGHz < 3 x 10 23 W Hz" 1 account 
for <10% of cluster heating; whereas the brightest two 
or three sources (including the BCG) contribute >90% 
of the total heating in proportion to their observed radio 
fluxes. 

However, a recent re-evaluation of the relationship be- 
tween jet power and r a dio luminos ity has been under- 
taken by IBirzan et all (|2004 l2008f) who find evidence 
that the magnetic field strength is not constant for these 



sources so that a shallower scaling is required: Pj e t cx 
L^- 35-0 - 7 . The largest uncertainty in these calculations 
appears to be due to source aging; the correlation for 
the youngest sourc es is at the steep end of the range 
( Birz anet al.l 12008) . Therefore, statistical results ob- 
tained from the RLF is more robust than a "snapshot" of 
a single cluster's AGN population at any given time. For 
example, while our modest-sized sample provides a mean 
radio galaxy census of ~ 2 per cluster at \j\aghz > 3 x 
10 23 W Hz" 1 , the spread in relative radio luminosities 
and thus ICM heating in these clusters is large. For ex- 
ample, 5 of these clusters have their total radio luminos- 
ity dominated by a BCG radio galaxy, 5 do not and one 
has no rad io galaxy above ou r flux limit at all. Neverthe- 
less, if the IBirzan et al.l (|2008f ) scaling is correct, weaker 
radio galaxies are proportionally more important than 
previously believed in heating the ICM. If Pj et oc L° 5 , 
sources weaker than 3 x 10 23 W Hz -1 account for more 
than half (~ 55%) of the heating and so must be taken 
into account in the modeling. This spatially-distributed 
(and moving!) heat input does however solve the prob- 
lem of heat distribution throughout the ICM. 

From the perspective of AGN heating models, the ra- 
dio galaxy data in Table [3] take on a new importance. 
Approximately half of these clusters have luminous ra- 
dio sources in their BCGs (the radio galaxy MS1008-R3 
is in a galaxy ~ 1 mag fainter than the BCG despite be- 
ing somewhat close to the ICM emission centroid) and 
one cluster (Abell 2111) has no luminous radio galaxies 
at all. If the AGN heating model for clusters is cor- 
rect, we might expect quite a wide variation in temper- 
ature profiles and/or temperature variations at different 
radii for clusters with and without luminous BCG radio 
sources. It is also possible that peculiar motions of these 
heat sources within clusters may smooth-out these dif- 
ferences. What is clear from the radio luminosities in 
Table [3] is that the luminous radio sources in non-BCGs 
cannot be ignored in cluster heating models. 

7. CONCLUSION 

We have presented a summary of first results from a 
large survey of X-ray and radio-selected AGN in clusters 
of galaxies. In the 11 clusters we chose as potential pro- 
genitors to the present-day Coma cluster, we find 20 ra- 
dio galaxies and 8 X-ray point sources (XPSs) with 75% 
of these AGN centrally concentrated within 500 kpc pro- 
jected radius. This central concentration of the AGN is 
significantly different (> 99.999% significance level) from 
the radial distribution of the bright (>L*) cluster red 
sequence (CRS) galaxies. This extreme central concen- 
tration strongly suggests that an AGN trigge ring mecha- 
nism s imilar to the Bondi accretion model of I Allen et alJ 
( 2006) is operating in these clusters. However, since most 
of these AGN are moving supersonically with respect to 
the cluster ICM, Hoyle-Littleton accretion must be op- 
erable, making any gas accretion onto a supermassive 
Black Hole extremely inefficient. Therefore, other mod- 
els for triggering and powering these radio-loud AGN 
should be considere d (e.g ., extracting Black Hole spin; 
iBlandford & Znai eld [1971 . 

Except for 3 X-ray /radio sources, cluster radio galax- 
ies and XPSs generally are not coincident within our ob- 
served luminosity limits of I^iaghz < 3 x 10 23 W Hz -1 
and L .3-8.ofceV > 10 42 ergs s _1 , respectively. While ra- 



10 



Hart, Stocke & Hallman 



dio galaxies are certainly associated with AGN activity, 
we questioned prima facie whether cluster XPSs in CRS 
ellipticals are truly AGN. These X-ray sources show no 
strong evidence for obscuration and are hosted by lumi- 
nous red galaxies with no evidence for typical Seyfert-like 
emission signatures. But non-AGN models for cluster 
XPSs (e.g., "cool core" emission in luminous galaxy nu- 
clei) seem quite implausible, with these XPSs being 1-2 
orders of magnitude brighter than expected in simula- 
tions and as obs erved in very nearby clusters like Virgo 
dSun et alJl200l . 

These XPSs are not spatially resolved in the Chandra 
images and have X-ray spectra which individually or in 
composite are better fit as power-law sources, not ther- 
mal sources. We conclude that the most viable model 
for the XPSs is that they are the beamed X-ray emis- 
sion from low-luminosity BL Lac Objects, similar to the 
hi gh-energy-pea ked B L Lacs (HBLs) studie d in detail 
bylStocke et all (119851) . iPadovani fe Giommil (|1995D , and 
iRector et al.l ([1999). HBLs are found in luminous, pas- 
sive ellipticals, have steep X-ray spectra with no evi- 
dence for internal absorption, optical spectra which in 
many cases are indistinguishable from normal giant el- 
lipticals and are found in rich groups and moderately 
rich clusters at similar reds hifts to those observed here 
(see lUrrv fc Padovanilll995l for a review of their proper- 
ties). The presence of excess blue light (blueward of the 
Ca II break) when compared to CRS elliptical galaxies 
in one XPS (Abell 963 X2) supports the HBL classifica- 
tion of these sources; a i?-band polarization detection in 
this object would solidify this assertion. One other XPSs 
with a smaller amount of excess blue light also is present 
in our sample. 

X-ray sources with similar properties have been termed 
X-ray emitting, passive elliptical g alaxies or "optically - 
dull AGN" in other contexts (e.g.. iBarger et al.ll200ibf ). 
and so we suggest based on our current observations that 
these sources are low luminosity HBLs as well. More 
accurate colors and calibrated spectra can be used to 
test this hypothesis by searching for the excess blue light 
we have found in at least one of our sample members. 
Therefore, both the XPSs and the radio sources are most 
readily identified as radio-loud AGN with jets which can 
transfer heat into the surrounding ICM. Extrapolations 
of the radio galaxy popul ation in these clusters using the 
observed RLF of AGN (jMachalski fc Godlowskil |2000T) 
and the XPS popula tion using the observed XLF of HBLs 
(|Rector et al.ll2000h converge to predict that 85-100% of 
all bright ellipticals within projected radii of 500 kpc 
from rich cluster centers will be detected by deep radio 
maps {h\AGHz > 10 21 - 4 W Hz -1 ) and X-ray imaging 
(detection limits of L .3_8.ofceV > 10 40 ergs s~ x ). There 
is already a hint th at this prediction is correct in that 
iSantra et al.1 (|2007| ) have found that all bright ellipti- 
cals in the very core of the Perseus Cluster have XPSs 
(L .3_8.ofceV = 10 38 ~ 40 ergs s _1 ) based on a very deep 
X-ray image of the core of that cluster. 

These results have an impact on models for ICM heat- 
ing. Based upon the RLF for AGN and the je t power vs. 
radio luminosity scaling law of lBicknelil (|1995f ). we expect 
that the few brightest radio galaxies (including the ra- 
dio source associated with the BCG) contribute >90% of 
the heating while the remainder at \j\aghz < 3xl0 23 W 



Hz -1 contribute <10%. But if the more recent scaling 
between jet power and radio luminosity advocated by 
iBirzan et alj (|2008| ) is correct, the lower power sources 
(below our X-ray and radio detection limits) contribute 
>55% of the heat input. Due to their relative motion, the 
non-BCG AGN heat the ICM more uniformly and would 
not be expected to possess obvious X-ray "bubbles" in 
the ICM associated with their jets due to their peculiar 
velocities. It seems likely that these ideas can be used to 
test whether AGN are respons ible for heating the clus - 
ter ICM as has been proposed (jRuszkowski et al.1 f2004). 
Putting together two of our results, (1) radio galaxies 
are more centrally-concentrated than CRS galaxies as a 
whole and (2) large numbers of low power radio galax- 
ies are responsible for most of the ICM heating, suggests 
that a feedback mechanism by which the density of the 
ICM triggers the heat sources to offset cooling may be 
operable in rich clusters of galaxies. We might also ex- 
pect that those clusters whose total radio luminosities 
are dominated by radio sources in the BCG might pos- 
sess a different temperature profile or detailed tempera- 
ture/density signatures than clusters dominated by non- 
BCG radio galaxies. The existence of these signatures 
and what these signa tures might be will b e addressed in 
a future publication (lHallman et al.ll2009l ). 
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FlG. 1. — Predicted X-ray tempe rature (keV) vers us redshift for Co ma Cluster progenitors. The model curves show results from detailed 
hydrodynamical simulations (similar to Evrard 1988 and Bialck ct al. 2001) of the evolution of clusters which evolve to masses similar to the Coma 
Cluster. Measures of Tx(z)/Tx{z-0) from the simulations are normalized by the Coma Cluster's present temperature to illustrate the variance 
of cluster properties. The solid line shows the median value for this distribution and the dashed line shows the 25th and 75th percentiles. Eleven 
clusters are chosen to be a representative sample of Coma Cluster progenitors at moderate rcdshifts (0.2<z<0.4). Sec § |3.1.1| for details on the bulk 
ICM temperature determination. 
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FlG. 2. — Composite color-magnitude diagram for eleven cluster fields. For each cluster, the color difference is computed relative to the average 
color of the cluster red sequence galaxies (see § |3.3.3t All optical sources within the 1 Mpc projected radius arc plotted with small dots. Brightest 
cluster galaxies (BCGs) are identified by an asterisk (with two BCGs for Abell 1758 which appears double-peaked in X-rays). Notice that cluster 
radio galaxies (filled triangles) and cluster X-ray point sources (filled circles) are located along the cluster red sequence, implying that they are 
hosted by passive, elliptical galaxies. Letter identifiers corresponds to X-ray point sources listed in Table [2] M*— -20.8 is marked by a vertical 
dashed line. The horizontal dotted lines arc representative of the 3<r spread (~ 0.36 mag) from the mean CRS color over the entire cluster sample. 
In general the cluster red sequence has a slope of ~ -0.03 based on SDSS (g.r) photometry of clusters in our sample, with more luminous CRS 
galaxies having redder colors. 
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FlG. 3. — The cumulative fraction of cluster red sequence (CRS) galaxies (dot-dash line), cluster radio galaxies (dotted line), cluster X-ray point 
sources (XPSs) (dashed line), and combined cluster radio galaxies and XPSs (solid line). The population of cluster radio galaxies and XPSs arc 
more centrally concentrated than cluster red sequence galaxies with ~ 75% of them located within a 500 kpc projected radius from the cluster X-ray 
centroid. A Kolmogorov-Smirnov test rules out at the 99.999% level that the radio galaxies and XPSs are drawn from the CRS parent population, 
suggesting that AGN triggering is connected to the ICM environment. 



XPSs and Radio Galaxies in Clusters of Galaxies 



15 



i — i — i — i — i — i — i — i — i — i — r 



C\J 
CO 

D) 
O 



CM 
CO 

O 



1.0 



o Cluster XPS 

■ Composite XPSs with 
radio non-detection 

A Composite XPSs with 
radio detection 



i — I 1 1 1 1 — I 1 1 1 1 — I — rr i , | — p — i 1 r 



0.5 



r=4.o 



# : 
O : 



r=3.o 



o.o 



r=2.o 



H 



r=i.o /. 



Or 



-0.5 



o Cr.v o 



<3 



j I j i i i I i i i i I i i : i -Jl I i i i i I i i i i I i i i l 



a 



51 = 0.5-0.9 keV 

52 = 0.9-2.5 keV 
H = 2.5-8.0 keV 



X 



-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 
C21 = log(S1/S2) 



FlG. 4. — Rest-frame X-ray colors of cluster X-ray point sources. The colors are defined as follows: SI— 0.5-0.9 keV, S2— 0.9-2.5 keV, and 
H=2. 5-8.0 keV. Overlaid is a grid of the expected color ratios in this color system for T— 1-4 and n^j— 0,0. 1,1.0, 10X 10 22 cm~ 2 . The open circles are 
our detected cluster XPSs. Letter identifiers correspond to the X-ray point sources listed in Table [2] The filled square is the composite color of 
the stacked XPS spectra without radio emission above our limits, excluding the one XPS (Abcll 963 XI, Object A) with a typical Scyfcrt optical 
emission line spectrum. The filled triangle is the composite color of the stacked XPS spectra with radio emission above our limits. The X-ray colors 
of the stacked XPS spectra are consistent with an unobscurcd AGN with T — 1.7. 
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FlG. 5. — X-ray (0.3- 8. OkeV) versus rad io (1.4 GHz) luminosity for cluster radio galaxies and X-ray point sources. Clust er sources are 
compared to typical FR Is (Donato ct al. 2004), ellipticals classified as "core" and "power-law" galaxies (Balmavcrdc & Capctti 2006), and Seyferts 
JCapetti iTBa lmavcrdc 2007). The vertical and horizontal dashed lines represent our radio power limit of T*xAGBz > 3xl0 23 W Hz 1 and X-ray 
luminosity limit (0.3-8.0 keV) of hx > 10 42 ergs s . Cluster radio galaxies are identified with filled triangles (radio-loud BCGs arc additionally 
identified with a black circle around the triangle), while cluster X-ray point sources arc identified with filled circles. Letter identifiers correspond to 
X-ray point sources listed in Tabic [2] Filled squares arc sources with both detectable X-ray and radio emission within our defined luminosity limits 
(with BCGs identified as above). Notice that our X-ray point sources with radio non-dctcctions and our radio galaxies with X-ray non-dctcctions 
do not overlap in this plot, suggesting two different populations of objects. 



TABLE 1 

"Road to Coma" Cluster Sample for z < 0.4 



Cluster 
Name 

fl \ 
W 


z 


Chandra 
ObsID 


ACIS 
chip 

W 


Exp. Time 
(ksec) 

/e:\ 

\ a ) 


T 

(keV) 

IP.\ 

W 


^ X , limit 

xl0~ 15 

(71 
(<) 


Lx, limit 

xlO 41 


VLA 
1.4 GHz Map 


Sl.4GHz,limit 

(mjy) 
(W) 


PlAGHz, limit 
XlO 23 


Comments 


MS 0440.5+0204 


0.197 


4196 


S3 


38.2 


6 1+ ' 6 

°- -0.5 


2.2 


2.1 


AS873 


0.1 


0.4 


AI72, A-array 


Abell 963 


0.206 


903 


S3 


35.8 


7 1+0.2 
'■ L -Q.2 


1.7 


2.1 


FIRST 


1.1 


1.4 




RX J0952. 8+5153 


0.214 


3195 


S3 


25.6 


c; 9+0.2 
°- z -0.2 


2.3 


3.1 


FIRST 


0.1 


0.5 




Abell 2111 


0.229 


544 


13 


9.7 


7 9+O.6 
—0 5 


8.2 


12.7* 


AS873 


0.3 


0.5 


FIRST 


MS 1455.0+2232 


0.260 


4192 


13 


75.6 


4 5+0-1 


2.0 


4.1 


AS873 


0.1 


0.8 


FIRST 


Abell 1758 


0.280 


2213 


S3 


35.8 


7 C-+0.5 
'•°-0.3 


1.9 


4.8 


FIRST 


0.4 


1.0 




MS1008. 1-1224 


0.306 


926 


13 


33.4 


6 1+ - 4 

°- -0.5 


3.2 


9.7 


Stocke et al. f 1999) 


0.6 


1.8 




MS2137.3-2353 


0.313 


5250 


S3 


23.5 




2.2 


7.0 


AB1022 


0.2 


0.6 




Abell 1995 


0.319 


906 


S3 


55.4 




1.2 


3.8 


AS873 


0.5 


1.6 


FIRST 


MS1358.4+6245 


0.330 


516 


S3 


34.1 


S.2t° .l 


1.5 


5.2 


FIRST 


0.5" 


1.7" 


Stocke et al. (1999) 


Abell 370 


0.373 


515 


S3 


59.3 




1.2 


5.8 


FIRST 


0.4 


1.9 


AL159, A-array 



Note. — Columns: (1) cluster name (2) cluster rcdshift (3) Chandra observation ID (4) ACIS aimpoint (5) Exposure time of Chandra observation after filtering for 
flaring events (6) ICM temperatures determined by XSPEC modeling of spectra extracted with a 1 Mpc (projected) radius and with cooling cores excised in some clusters. 
Sec § 13. 1. 11 for details. (7) Flux limit (0.3-8.0 keV) in ergs cm" 2 s" 1 for a point source near the edge of our survey region (R — 1 Mpc). Sec § |3.1.2| for more details. (8) 
Luminosity limit (0.3-8.0 keV) in ergs s — 1 for the flux limits in Column 7 (9) VLA program or Reference; AS873 Observations were obtained in B or C-array configuration 
(10-11) 1.4 GHz Radio flux density in mjy and power limits in W Hz -1 for a 3er detection at the edge of our survey region (12) Additional VLA programs used to identify 
radio sources 

The X-ray luminosity limit in Abell 2111 is slightly greater than our minimum threshold of Lx > 1 X 10 42 ergs s _1 . Sec § |3.1.2| for our justification on its inclusion in 
the cluster sample. These limits are extrapolated from 5 GHz observations assuming F, y cx v~ a and a— 0.7. Sec § |3.2| for details. 



TABLE 2 

Cluster X-ray Point Sources with L 3 _ s . 0keV > 10 42 ergs s _1 



Object 


a 


5 


2 


r 


Radius 


Net Counts 




^X,rest x 10 42 


Comments 




(J2000) 


(J2000) 




mag 


(hr„ Mpc) 
V TO 


(0.3-8.0 keV) 


(ergs cm - " s~ 1 ) 


ferss s — ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(A) A963-X1 


10:17:00.8 


+39:04:33.0 


0.209 a 


18.22 


0.360 


2150.0 ± 48.8 


394.2 ± 8.9 


45.8 ± 1.0 


[O III], Ha emission lines 


(B) MS1008-X3 


10:10:35.3 


-12:40:22.1 


0.3095 6 


19.39 


0.264 


60.2 ± 9.5 


23.4 ± 3.7 


6.5 ± 1.0 




(C) A370-X6 


02:40:00.2 


-01:32:34.1 


0.380 ac 


19.94 


0.831 


62.3 ± 9.5 


6.5 ± 1.0 


2.9 ± 0.4 




(D) A963-X2 


10:16:58.7 


+39:02:12.7 


0.206" 


20.74 


0.237 


115.1 ± 12.8 


23.8 ± 2.7 


2.8 ± 0.3 


Passive spectrum with blue excess 


(E) MS1455-X2 


14:56:58.7 


+22:18:46.3 


0.2582 6 


20.05 


1.008 


26.0 ± 6.5 


7.7 ± 1.9 


1.5 ± 0.4 


Coincident with radio source, MS1455-R1; 


















Probable blue excess 


(F) A1758-BCG 


13:32:38.3 


+50:33:36.2 


0.279" 


17.19 


0.000 


30.9 ± 9.9 


5.7 ± 1.8 


1.3 ± 0.4 


Coincident with radio source, A1758-R4 


(G) MS0440-X8 


04:43:14.2 


+02:12:04.7 


0.1996 6 


18.37* 


0.403 


48.6 ± 8.4 


11.0 ± 1.9 


1.2 ± 0.2 


Coincident with radio source, MS0440-R5 


(H) RXJ0952-X6 


09:52:36.2 


+51:56:32.9 


0.22° 


19.47 


0.842 


33.4 ± 7.8 


8.7 ± 1.9 


1.1 ± 0.2 





Note. — Columns: (1) Cluster X-ray point source IDs, with same letter ID as displayed in Figs. \2\ & \E\ (2-3) (J2000) RA 8z DEC of X-ray point source (4) redshift (5) observed Sloan 
r-band magnitude, except for MS0440-X8 which is Gunn r-band magnitude (6) projected radial distance (Mpc) from the peak in the cluster X-ray emission (7) net counts in 0.3-8 keV bandpass 
(8) 0.3-8.0 keV Flux (9) 0.3-8.0 keV Rest-frame Luminosity assuming a power-law spectrum with a photon index of 1.7 (N_e oc E — r ). 

Gunn r-magnitudc a This work b CNOC; [Yee et all 1 1998ft c [Barger et all d2001a|) 



TABLE 3 

Cluster Radio Galaxies with Piaghz > 3 x 10 23 W Hz -1 









c 



z 


r 


(g- r ) 


(r-i) 




CRS 






Radius 




riAGHz X 


i n23 
1U 


L x xl0 42 




U bj ect 


(J2000) 




(J2000) 




mag 


color 


color 




color 






Mpc 


(mJy) 


(WT IT.. — 

(W Hz 


) 


ergs s 


i 




(!) 


I 2 ) 






(4) 


(.5) 




(') 




(8) 






(9) 


(10) 






(12) 






MS0440-R3* 


04:43:10 


.0 


+02:10:18.4 


0.199 a 


17.03** 


0.99 




1 


.00 


± 





.30 


0. 


.001 


30.2 ± 0.07 


34.2 ± 


0.1 






4 


MS0440-R4 


04:43:09 


.4 


+02:10:00.0 


0.193° 


18.67** 


0.93 




1 


.00 


± 





.30 





.062 


7.5 ± 0.08 


8.0 ± 


0.1 




0. 


2 


MS0440-R5 


04:43:14 


.2 


+02:12:04.2 


0.120° 


18.37** 


0.92 




1 


.00 


± 





.30 





.392 


2.7 ± 0.2 


3.1 ± 


0.2 


1.2 ± 


0. 


2 


A963-R2 


10:17:22 


.0 


+39:04:39.1 


0.197 c 


17.14 


1.34 




1 


.31 


± 


0. 


.21 


0. 


.806 


3.82 ± 0.38 


4.2 ± 


0.4 


< 


0. 


2 


RXJ0952-R1* 


09:52:49 


.1 


+51:53:05.0 


0.215 d 


16.56 


1.28 




1 


.31 


± 


0. 


.18 


0. 


.009 


13.70 ± 0.13 


18.5 ± 


0.2 




1. 


5 


RXJ0952-R3 


09:52:46 


.9 


+51:54:42.2 


0.22 c 


19.48 


1.44 




1 


.31 


± 





.18 





.351 


5.54 ± 0.13 


7.9 ± 


0.2 


< 


0. 


1 


MS1455-R1 


14:56:58 


.7 


+22:18:46.3 


0.258-? 


20.05 


1.32 




1 


.47 


± 





.30 


1 


.014 


19.5 ± 0.1 


39.9 ± 


0.2 


1.5 ± 


0. 


.4 


MS1455-R2* 


14:57:15 


.1 


+22:20:34.3 


0.258-? 


16.50 


1.24 




1 


.47 


± 





.30 


0. 


.007 


15.9 ± 0.1 


32.5 ± 


0.2 


< 


2. 


3 


MS1455-R3 


14:57:08 


.5 


+22:20:14.5 


0.25 ° 


18.39 


1.56 




1 


.47 


± 





.30 


0. 


.379 


14.2 ± 0.1 


27.0 ± 


0.2 


< 


0. 


2 


A1758-R4* 


13:32:38 


.4 


+50:33:35.9 


0.278 c 


17.19 


1.55 




1 


.52 


± 





.39 





.000 


6.03 ± 0.15 


14.6 ± 


0.1 


1.3 ± 


0. 


1 


A1758-R6 


13:33:02 


.0 


+50:29:28.6 


0.279 d 


17.96 


1.60 




1 


.52 


± 





.39 





.670 


1.76 ± 0.15 


4.3 ± 


0.1 


< 


0. 


9 


MS1008-R1 


10:10:31 


.1 


-12:37:02.6 


0.311°° 


18.26** 


1.19 




1 


.25 


± 





.40 


0. 


.730 


25.6 ± 0.2 


80.3 ± 


0.6 


< 


0. 


3 


MS1008-R2 


10:10:29 


.4 


-12:38:27.4 


0.300 c 


18.72** 


1.23 




1 


.25 


± 





.40 


0. 


.393 


12.9 ± 0.2 


37.2 ± 


0.6 


< 


1. 


1 


MS1008-R3 


10:10:32 


.3 


-12:39:34.1 


0.301^ 


18.75** 


1.15 




1 


.25 


± 





.10 





.042 


5.9 ± 0.2 


17.2 ± 


0.6 


< 


3. 


1 


MS2137-R4* 


21:40:15 


.1 


-23:39:40.0 


0.314° 


17.36 




0.57 





.50 


± 





.25 


0. 


.001 


3.8 ± 0.5 


12.2 ± 


1.6 


< 


8. 


.3 


A1995-R1 


14:53:00 


.5 


+58:03:20.0 


0.325° 


19.15 




0.58 





.59 


± 





.21 





.151 


8.43 ± 0.15 


29.3 ± 


0.5 


< 


0. 


.3 


MS1358-R1* 


13:59:50 


.6 


+62:31:04.2 


0.326° 


17.75 




0.66 





.57 


± 





33 


0. 


.005 


2.1 ± 0.9 


7.4 ± 


3.2 


< 


3. 


6 


A370-R1 


02:39:55 


.3 


-01:34:05.4 


0.382°'' 


19.97 




0.68 





.75 


± 





.33 





.224 


4.6 ± 0.15 


23.3 ± 


0.8 


< 


0. 


.4 


A370-R2 


02:39:52 


.9 


-01:35:01.4 


0.376° 


21.21 




0.88 





.75 


± 





.33 





.123 


1.51 ± 0.15 


7.4 ± 


0.7 


< 


0. 


3 


A370-R3 


02:39:56 


.3 


-01:34:29.3 


0.370' 1 


19.68 




0.65 





.75 


± 





.33 


0. 


.280 


1.36 ± 0.05 


6.4 ± 


0.2 


< 


0. 


3 



Note. — Columns: (1) Radio Galaxy ID with brightest cluster galaxies (BCGs) denoted by an asterisk. For Abell 2111, no radio galaxy was detected 
above our radio power limit. (2-3) J2000 Coordinates (4) object redshifts (see reference below) (5) apparent Sloan r-band magnitude except for those with 
double asterisks which are Gunn r-band magnitudes (6) Observed color (g-r) or (7) (r-i) (8) Mean color of the cluster red sequence. The mean Sloan color 
and 3er spread of the cluster red sequence for M* < —20.8. (9) Projected radial distance from the X-ray peak emission. Note that the X-ray image of 
Abell 1758 is a double-peaked in appearance, indicative of a merging system. Therefore, the projected radii values arc relative to the NW and SE X-ray 
emission clumps. (10) observed flux density of the source (11) observed radio power of the source (12) Rest frame X-ray Luminosity (0.3-8 keV) assuming 
a power-law spectrum with T— 1.7 (N e oc E~ r ). The BCG X-ray flux limits are calculated assuming a minimum of 3cr above the diffuse cluster emission. 
For the remaining radio sources, X-ray flux limits we re determin ed using ne t coun ts (observed minus local backgro und) plus 3(7 above the bac kground noise. 
Refere nce s for source redshif ts: C alGioia et al.Hl998fl (b) CNOC: lYee et al. | fl998T) (c) This work (d) SDSS DR6 (e) IShectman et all {19961) Cfl lMartini et al.l 
(120071) fgl iPatel eraD l2000h ("hi lMellier et al.1 119881) fii lStocke et al.l dl999h (j) CNOC; communication from E. Ellingson 

Radio source is hosted in the Brightest Cluster Galaxy Gunn r-band magnitude 
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TABLE 4 

AGN Fraction for Cluster Red Sequence (CRS) Galaxies 



Cluster 


^* c as 


N r 


Ir 


N x 


fx 


Color Used 


Survey 


(1) 


(2) 


(3) 


(J) 


(5) 




(7) 




(8) 


MS 0440.5+0204 


18 


3 a 


0.167 


1 


0.056 


Gunn 


(g-r) 


CNOC 


Abell 963 


18 


1 


0.021 


2 


0.042 


Sloan 


(g-r) 


SDSS DR6 


RX J0952.8+5153 


23 


2 a 


0.087 


1 


0.044 


Sloan 


(g-r) 


SDSS DR6 


Abell 2111 


72 





0.000 





0.000 


Sloan 


(g-r) 


SDSS DR6 


MS 1455.0+2232 


46 


3 a 


0.065 


1 


0.022 


Gunn 


(g-r) 


CNOC 


Abell 1758 


84 


2 a 


0.024 


1 


0.012 


Sloan 


(g-r) 


SDSS DR6 


MS1008. 1-1224 


69 


3 


0.044 


1 


0.015 


Gunn 


(g-r) 


CNOC 


MS2137.3-2353 


56 


1" 


0.018 





0.000 


Sloan 


'r-i) 


ChaMP 


Abell 1995 


71 


1 


0.014 





0.000 


Sloan 


'x-ij 


SDSS DR6 


MS1358. 4+6245 


67 


l a 


0.015 





0.000 


Sloan 


'r-i) 


SDSS DR6 


Abell 370 


111 


3 


0.027 


1 


0.009 


Sloan 


'r-i) 


APO 2007 


Average 






0.044 




0.018 








Sum 


665 


20 


0.030 


8 


0.012 









Note. — Columns: (1) Cluster name (2) number of CRS galaxies with M* < —20.8 and 
projected radius < 1 Mpc from cluster center as estimated using the procedure described 
in § 13.3.31 (3) number of radio galaxies with Af* < -20.8 and Pi.iGXz > 3 x 10 23 W 
Hz -1 (4) radio galaxy fraction within our limits (5) number of X-ray point sources with 
M* < —20.8 and Lq .3_gfceV > 10 42 ergs s _1 (6) X-ray point source fraction within our 
limits (7) observed color used for the CRS galaxy estimation (8) source of the imaging 
data 

a The brightest cluster galaxy (BCG) is a radio source. 



TABLE 5 

Comparison of Observed vs. Expected Optical Colors for 2 "Blue" XPSs 



Object 


Alternative ID 


Color 


Observed 


CRS a 


Comments 


Abell 963 X2 


Object D 


(u-g) 


0.70 ± 0.13 


1.9 ± 0.5 


APO 3.5m 






(g-r) 


0.77 ± 0.06 


1.26 ± 0.07 


APO 3.5m 


MS1455 X2 


Object E 


( u -g) 


> 0.6 


2.2 ± 1.1 


SDSS DR6 






(g-r) 


1.32 ± 0.08 


1.45 ± 0.04 


SDSS DR6 



a The expected mean color (SDSS DR6) of a cluster red sequence (CRS) galaxy of 
similar luminosity to the object and its associated lcr spread about the mean. 



